Molecular Line Observations of a Carbon- Chain- Rich Core L492 
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C<l : ABSTRACT 

vn . 

O ! We report on molecular abundances and distributions in a starless dense 

Q I core L492. We have found that the abundances of carbon-chain molecules such 

"^ I as CCS, C3S, HC3N, HC5N, and HC7N are comparable to those in chemically 

Ph| young dark cloud cores called "carbon-chain-producing regions" , such as L1495B, 

O ' L1521B, L1521E, and TMC-1. This is the first dark cloud core with extremely 

'^ I rich in carbon-chain-molecules that is found outside the Taurus region. In ad- 

f^ I dition, the deuterium fractionation ratios of DNC/HNC and DCO"'"/HCO"'" are 

.^ ' also comparable to those in carbon-chain-producing regions, being significantly 

/\f ' lower than those in the evolved prestellar cores such as L1498 and L1544. On 

c^ ' the other hand, the abundances of NH3 and N2H"'" are systematically higher than 

those in carbon-chain-producing regions. Our mapping observations reveal that 
the central hole of molecular distributions, which were reported for CCS and C^^S 
in evolved prestellar cores is not significant in L492, indicating that the depletion 
factor of molecules is not very high. Furthermore, L492 is dynamically more 
evolved than carbon-chain-producing regions, and the protostellar collapse has 
started like L1498 and L1544. Therefore, it is likely that the chemical and dynam- 
ical evolutionary stage of L492 is intermediate between carbon-chain-producing 
regions (L1495B, L1521B, L1521E, and TMC-1) and evolved prestellar cores 
(L1498 and L1544). 

Subject headings: ISM: abundances — ISM:individual(L492) — ISM:Molecules — 
radio lines: ISM 
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Introduction 



In the last decade, special attention has been paid to starless cores or prestellar cores 
in dark clouds in order to investigate initial conditions of protostellar collapse. Since the 
detection of infalling motion in a prestellar core L1544 in the Taurus Molecular Cloud (Myers 
et al. 1996), a number of observational studies on L1544 and similar prestellar cores such as 
L1498 in the Taurus Molecular Cloud have been carried out with molecular lines and dust 
continuum emission (e.g. Tafalla et al. 2002, 2004). One of the most interesting features 
found in L1498 and L1544 is an evidence of molecular depletion in the central part of the 
cores (e.g. Kuiper et al. 1996; Ohashi et al. 1999; Caselli et al. 2002b). Both observational 
and theoretical studies reveal that abundances of CO, CS, and CCS are depleted in the 
chemically evolved cores, while NH3 and N2H+ are not (Bergin & Langer 1997; Aikawa et 
al. 2005). These studies have promoted understanding of chemical properties as well as 
kinematics of evolved prestellar cores just under the collapsing phase. 

Although detailed studies on evolved prestellar cores were extensively carried out, no 
systematic studies on younger cores have been made for a long time, except for TMC-1 (e.g. 
Hirahara et al. 1992, 1995; Pratap et al. 1997). It is because most of dense cores observed in 
the previous studies were based on the catalog of Benson & Myers (1989), which are prepared 
by the NII3 observations. On the other hand, Suzuki et al. (1992) pointed out that NII3 is not 
always a good tracer of dense cores because of the chemical abundance variation from core 
to core. They found anticorrelation between abundances of NH3 and carbon-chain-molecules 
such as CCS, HC3N, and HC5N, and identified a few starless cores where the carbon-chain 
molecules are abundant while NH3 is deficient. They are L1495B, L1521B, L1521E, and the 
cyanopolyyne peak of TMC-1, and are called "carbon-chain-producing regions". Chemical 
model calculations suggest that the carbon-chain-producing regions are in the early stage of 
chemical evolution while dense cores traced by the NH3 lines are rather evolved (Suzuki et 
al. 1992; Bergin & Langer 1997). 

In order to investigate basic physical and chemical properties of carbon-chain-producing 
regions, we carried out detailed observations of L1521E (Hirota et al. 2002), L1495B and 
L1521B (Hirota et al. 2004). For each of these three regions, we found a compact dense 
core traced by the H^^CO^, CCS, C3S, and HC3N lines and their distributions have a single 
peak at the same position. Such a distribution is different from those in evolved prestellar 
cores, L1498 and L1544. In addition, abundances of NH3 and N2H"'" are extremely lower in 
these three regions than in other prestellar cores, while those of carbon-chain-molecules are 
systematically higher. According to our survey of deuterated molecules, deuterium fraction- 
ation ratios of DNC/HNC and DCO+/HCO+ are systematically lower in L1495B, L1521B, 
and L1521E than in the other dark cores (Hirota et al. 2001, 2003). These features suggest 
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that L1495B, L1521B, and L1521E would be in the early stage of dynamical and chemical 
evolution, and the degree of depletion of molecules is possibly lower than in other evolved 
cores (Tafalla & Santiago 2004). 

Although detailed studies on carbon-chain-producing regions are important for under- 
standing of chemical and physical evolution of dense cores, only four carbon-chain-producing 
regions were identified in Taurus Molecular Cloud. In order to search for carbon-chain- 
producing regions outside the Taurus Molecular Cloud, we carried out a survey of CCS, 
HC3N, and HC5N toward 31 dark cloud cores. As a result, we detected a possible candidate 
for another carbon-chain-producing region L492 in the Aquila Rift (Hirota et al. in prepara- 
tion). L492 is one of the prestellar cores identified as a strong infalling candidate (Lee et al. 
2001), and it is observed with the dust continuum emission as well as the N2D"'" and N2H"'" 
lines (Crapsi et al. 2005). In this paper, we report on the results of extensive molecular line 
observations of L492. 



2. Observations 

We carried out molecular line observations of L492 with the 45 m radio telescope at 
Nobeyama Radio Observatory (NRO)^ in February 2003 and March 2004. The observed 
lines are summarized in Table 1. 

We used cooled HEMT receivers for the lines in the 22-23 GHz band, SIS mixer receivers 
for the lines in the 45-48 GHz and 72-96 GHz bands, the SIS 25-beam array receiver system 
(BEARS) for the ^^CO and C^^O lines in the 109-110 GHz band. The main-beam efficiencies 
{r]^b) were 0.8, 0.7, 0.5, and 0.5 for the 22-23 GHz, 45-49 GHz, 72-96 GHz, and 109-110 GHz 
bands, respectively. The beam sizes were 73", 36", 20", and 16" for the 22-23 GHz, 45-49 
GHz, 72-96 GHz, and 109-110 GHz bands, respectively. Acousto-optical radio spectrometers 
with the frequency resolution of 37 kHz were used for the backend except for the ^^CO and 
C^^O observations, in which autocorrelators with the frequency resolution of 31.25 kHz were 
used. 

We took the reference position of L492 to be 02000 = 18^*15^46^1, ^2000 = -03°46'13"(Lee 
et al. 2001; Crapsi et al. 2005). Observations in the 22-23 GHz and 45-49 GHz bands were 
performed in the position-switching mode. The off position was taken to be 10' away from 
the source position, which is outside the core of L492 with the radius of 3.2'(Lee & Myers 



^Nobeyama Radio Observatory is a branch of the National Astronomical Observatory of Japan, an in- 
tcruniversity research institute operated by the Ministry of Education, Science, Sports and Culture of Japan 
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1999). Observations in the 72-96 GHz and 109-110 GHz bands were performed in the fre- 
quency switching mode with the frequency offset was set to be 7.685 MHz. The grid spacing 
in the mapping observations was taken to be 40-80", 40", and 20" for the 22-23 GHz, 45-49 
GHz, and 72-96 GHz bands, respectively, as shown in Figure 1. Since the beam separation of 
the BEARS is 41.1", the grid spacing of the ^^CO and C^^O maps were taken to be 20.55". 
Note that we did not carry out the Niquist samphng because of hmited observing time. 
Pointing was checked by observing a nearby SiO maser source, R-Aql, every 1-2 hours, and 
the pointing accuracy was estimated to be better than 5". The antenna temperature was 
cahbrated by a usual chopper-wheel method. 



3. Results 

3.1. Molecular Distributions and Velocity Structure 

Figure 1 shows the integrated intensity maps, and Figures 2-5 show an example of the 
observed spectrum for each molecule. We only show the maps of the lines marked in Table 
1. The molecular gas is distributed along the north-south ridge traced by the ^^CO and 
C^^O lines, and the dense core of L492 traced by other lines is located at the southern edge 
of the ridge. Most of the maps show quite similar structure with a single peak and a simple 
elliptical shape with the average diameter of 100", corresponding to 0.1 pc at the distance 
of 200 pc (Lee & Myers 1999). Although the peak positions are different from each other, 
this is partly due to the different beam size and grid spacing. The single peak distribution 
of the C^^S and CCS in L492 clearly contrasts with the double peak distribution in L1498 
and L1544 (Kuiper et al. 1996; Ohashi et al. 1999; Tafalla et al. 2002, 2004), and is rather 
similar to the distribution in the carbon-chain-producing regions such as L1495B, L1521B, 
and L1521E (Hirota et al. 2002, 2004). The peak position of the N2H+ map is 20" south 
of those of the dust continuum, N2D"'", and N2H"'" maps reported by Crapsi et al. (2005), 
probably due to insufficient spatial resolution or coarse grid spacing. 

On the other hand, we found a marginal double peak structure in the maps of the J=l-0 
line of H^^CO"*" and the 75-65 line of CCS, although the northern peak has the statistical 
significance of only 2 a level. Because the optical depth of the J=l-0 line of H^^CO^ and 
the 73-65 line of CCS are 0.82 and 0.54, respectively, as discussed later, they would not be 
suffered from self- absorption. A possible explanation is depletion of H^^CO"'" and CCS in the 
central part of the core. Note that the double peak structure is not evident in the CCS 43-82 
map as mentioned above. Since the critical density for the 70-65 liiie is higher than that of 
the 43-82 line, we could detect the depletion of CCS only with the 70-65 line. If the double 
peak structure is real, the depletion of molecules would have just started at the central part 
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of L492, which can be seen only in the 73-65 hne of CCS and the J=l-0 hne of H^^CO^. 
Further high resolution observations would be necessary to confirm this. 

It is interesting to compare the deuterium fractionation ratio within the dense core, 
because the deuterium fractionation ratio refiects chemical evolutionary stage and degree 
of depletion (Caselli et al. 2002b; Hirota et al. 2001, 2003). Figure 6 shows the integrated 
intensity ratios of DCO^/H^^CO^ and DNC/HN^^C as a function of the distance from 
the dust continuum peak. Because these lines are optically thin as discussed later, the 
integrated intensity ratios approximately represent the abundance ratios. The maximum 
DCO'''/H^^CO^ and DNC/HN^^C ratios are 1.6 and 2.2, respectively, and their positions 
are offset from the DCO^ and DNC peaks and also from the dust continuum peak. The 
deuterium fractionation ratios are not significantly enhanced in the central part of the core. 
We might be able to see an only slight increase in the DCO^/H^^CO^ ratio within the inner 
40" area. However, this is not so clear as in the case of L1544 (Caselli et al. 2002b). As for 
the DNC/HN^^C ratio, such an enhancement is not seen, being consistent with other cores 
(Hirota et al. 2003). 

We also investigate velocity structure of the core using the CCS (70-65) line, because 
this line has no hyperfine splitting and was observed with a relatively high velocity resolution 
(0.136 km s~^). Figures 7 and 8 show the channel maps and position-velocity diagrams of 
the CCS (73-65) line, respectively. We found a global velocity gradient of 3 km s~^ pc~^ 
in the east-west direction across the core. On the other hand, we could not find other 
significant change in the linewidth (e.g. Crapsi et al. 2005). According to Lee et al. (2001), 
the optically thick line, CS(J=2-1), toward L492 has the asymmetric line profile due to an 
infalling motion. However, either double peak or asymmetric line profiles are not detected 
in our observations because we observed only optically thin lines. 



3.2. Abundances for Individual Molecules 

In order to evaluate abundances of the observed molecules, line parameters for all the 
observed lines were determined by the Gaussian fit, as summarized in Table 1. We ignored 
unresolved hyperfine structures of the NH3 and HC3N lines, and hence, their linewidths are 
shown to be broader than the others. The parameters at the (0", 0") position for NH3 and 
HC7N (22-23 GHz band), the (40",0") position for other carbon-chain molecules (22-23 GHz, 
45-49 GHz and 81-96 GHz bands), and the (20", 20") position for the rest (72-93 GHz and 
109-110 GHz bands) are listed in Table 1. These positions correspond to the peak positions 
of the integrated intensity map for each molecule. Using these line parameters, we calculated 
the column densities of the observed molecules by the consistent way employed in previous 
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works (e.g. Suzuki et al. 1992; Hirota et al. 1998, 2001, 2002, 2004) in order to compare 
the present results with those for L1495B, L1521B, L1521E, and TMC-1. The methods are 
briefly described below, and the results are summarized in Tables 2 and 3. The references for 
collision rates and dipole moments are also summarized in Tables 2 and 3. In the analysis, 
the kinetic temperature in L492 is assumed to be 10 K, which was derived from the NH3 lines 
in the present study. Because we observed molecular lines with the different grid spacing 
depending on the beam size, we can not compare the line parameters and column densities 
of all the observed molecules at the same position. Therefore, we list the line parameters and 
molecular column densities at the above positions. Although the peak positions are different 
from each other, this is partly due to the different beam size and grid spacing. Therefore, 
the column densities for the positions listed in Tables 2 and 3 would be those for the same 
volume of gas. According to our maps shown in Figure 1, the difference in the integrated 
intensities due to the position difference is less than 20%, and hence, the uncertainty in the 
column densities due to the position difference is estimated to be 20%. 



3.2.1. C^^S 

We derived the column density of C^^S and the H2 density using the large velocity 
gradient (LVG) model (Goldreich & Kwan 1974). The method is the same as that adopted 
in Hirota et al. (1998). We ignored the difference in the beam sizes between the J=l-0 and 
2-1 lines in the analysis. Because the core size of L492 (100") is sufficiently larger than the 
beam sizes for the J=l-0 (36") and J=2-l line (20"), this approximation does not seriously 
affect the results (e.g. Hirota et al. 1998). The H2 density is derived to be 9.0x10"^ cm"'^. 



3.2.2. CCS 

We derived the column density of CCS and the H2 density using the LVG model. Because 
the beam size of the JAr=2i-lo lines are as large as the core size, beam dilution effect would 
be significant, and hence, we only used the results of the JAr=43-32 and JAr=76-65 lines. For 
these lines, we ignored the difference in the beam sizes of the two transitions, as in the case 
of the C'^^S analysis. The H2 density is derived to be 6.8x10'^ cm"^. 

For comparison, we carried out LTE analysis assuming the excitation temperature of 5 
K (Suzuki et al. 1992; Hirota et al. 1998, 2002, 2004). In this case, column density of CCS is 
three times larger than that derived from the LVG analysis. Therefore, the column densities 
of CCS for other cores obtained by the LTE analysis, in which the excitation temperature 
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is fixed to 5 K, would be overestimated. In this paper, the LTE value is employed for 
consistency with previous results. 



3.2.3. DCO+, DNC, H^^'CO+, HN^^C 

We derived the column densities of DCO+, DNC, YL^^GO+ , and HNi^C by the LVG 
calculations. The method is described in Hirota et al. (2001). Because we did not carry out 
multi-transition observations for these molecules, we assumed the H2 density to be 9.0x10*^ 
cm~^, which is derived from the C^^S data. If we assume the H2 density to be 2.0x10^ cm~^ 
(Crapsi et al. 2005), the column densities of DCO"^ and H^^CO"'" decrease by 10-25%, while 
those of DNC and HN^^q decrease by a factor of 2. In both case, the DC0+/Hi3C0+ and 
DNC/HN^^C ratios change by less than 10%. 



3.2.4. ^^co, a^o 

The column densities of ^^CO and C^^O were derived by assuming the ^^CO/C^^O ratio 
of the local ISM value, 7.3 (Wilson & Rood 1994), and the common excitation temperature. 
The excitation temperature was determined to be 8.3 K. 



3.2.5. HC3N 

At first, we carried out the LVG analysis using the J =5-4 and 9-8 lines to derive the 
column density of IIC3N and the H2 density. However, the fitting did not converge prob- 
ably because of the large optical depth of the J=5-4 line. Because the satellite hyperfine 
component of the J=5-4 line is clearly detected as shown in Figure 3, we can determine the 
optical depth and the excitation temperature of the J=5-4 line of HC3N. The total optical 
depth and the excitation temperature were calculated to be 8.8 and 6.4 K, respectively. The 
excitation temperature determined here is consistent with that employed by Suzuki et al. 
(1992), 6.5 K. Then we determined the column density of HC3N from the total optical depth 
assuming the LTE condition with the excitation temperature of 6.4 K. 



3.2.6. CsS, HC5N, HC7N 

For C3S, HC5N, and HC7N, we determined their column densities by the consistent way 
adopted in the previous hteratures (e.g. Hirota et ah 2002, 2004; Cernicharo et al. 1986). 
We assumed that the excitation temperatures of 5.5 K for C3S, 6.5 K for HC5N, and 10. OK 
for HC7N. 



3.2.7. N2H+ 

Since all the hyperfine components of the N2H+ line are detected as shown in Figure 4, 
we were able to determine the optical depth and the excitation temperature of N2H"'" to be 
6.1 and 4.9 K, respectively. The excitation temperature determined here is consistent with 
those reported by Caselli et al. (2002a), 5 K. Then we derived the column density of N2H"'" 
assuming the LTE condition with the excitation temperature of 4.9 K. 



3.2.8. NH3 

Using the intensity ratios of hyperfine components of the (1,1) lines, we derived the 
excitation temperature to be 4.1 K. In addition, we determined the rotation temperature 
of the para NH3 to be 9.5^-^' K by the methods described in Ho & Townes (1983). Thus, 
we calculated the column density of NH3 assuming the LTE condition with the rotation 
temperature of 10 K and the ortho-to-para ratio of 1. 



4. Discussions 

The column densities of the observed molecules are summarized in Table 4. For com- 
parison, the corresponding column densities in the carbon-chain-producing regions are also 
listed in Table 4. The column densities of the carbon-chain molecules in L492 are significantly 
higher than those in typical dark cloud cores (Suzuki et al. 1992) and are comparable to those 
in the carbon-chain-producing regions such as L1495B, L1521B, L1521E, and TMC-1. Since 
the carbon- chain-producing regions have so far been recognized in the Taurus Molecular 
Cloud, L492 is the first dark cloud core with extremely high abundances of carbon-chain- 
molecules that is found outside the Taurus Molecular Cloud. On the other hand, the column 
densities of NH3 and N2H+ in L492 are higher than those in the carbon-chain-producing re- 
gions, and are comparable to those in typical dark cloud cores (Benson & Myers 1989; Suzuki 



-9- 



et al. 1992; Caselli et al. 2002a). We stress that the column densities in L492 are almost 
comparable to those of TMC-1 for most of the observed molecules, except for the slightly 
high column density of NH3 and the low column densities of sulfur-bearing carbon-chain 
molecules and HC7N. Thus the chemical composition in L492 is close to that of TMC-1. 

In order to discuss chemical evolutionary stages of L492 and other starless cores, we 
summarize the selected molecular abundance ratios which are thought to be used as a chem- 
ical clock, as shown in Table 5. The most popular indicator is the NH3/CCS ratio (Suzuki et 
al. 1992). The NH3/CCS ratio of 6.5 in L492 is intermediate between those of carbon-chain- 
producing regions (2.6-3.8) and evolved prestellar core such as L1498 and L1544 (15-25). For 
consistency, we employed the NH3/CCS ratios in L1498 and L1544 reported by Suzuki et al. 
(1992). Uncertainties in the NH3/CCS ratios are quite large due to difference in telescope 
beam sizes, observed positions, and methods of data analysis. In fact, different values are 
reported by Aikawa et al. (2005) (14-25 and 9.0 for L1498 and L1544, respectively). If we 
use the result of the LVG calculations for the column density of CCS, the NH3/CCS ratio 
in L492 is 21, which is larger by a factor of 3 than that of the LTE value. Even such uncer- 
tainties are considered, we can say that the NH3/CCS ratio in L492 is intermediate or close 
to those of L1498 and L1544 rather than those in carbon-chain-producing regions. 

Another good indicator for chemical evolutionary stage is a deuterium fractionation 
ratio. Because deuterium fractionation proceeds through the exothermic isotope exchange 
reaction of B.^ and HD in the gas-phase, the 1120"*" /H^ ratio gradually increase with time 
(e.g. Roberts & Millar 2000; Turner 2001; Hirota et al. 2001; Saito et al. 2002). In addition, 
depletion of molecules make the life time of II2D+ longer, resulting in a high deuterium 
fractionation ratio of H2D+/H^. Consequently, a deuterium fractionation ratio is thought 
to be enhanced in evolved cores (e.g. Roberts & Millar 2000; Turner 2001; Hirota et al. 
2001; Saito et al. 2002). The deuterium fractionation ratio, DNC/HN^^C, of 1.27 in L492 
is comparable to those of the carbon-chain-producing regions (0.66-1.25), and is smaller 
than those in L1498 (1.91) and L1544 (3.0). The DC0+/Hi3C0+ ratio of 0.80 is also 
comparable to those of the carbon-chain-producing regions (0.63-1.10) and smaller than 
those in L1498 (2.7) and L1544 (3.1-9.2). The lower DC0+/Hi3C0+ ratio in L492 suggests 
that molecular depletion is not very significant in comparison with L1544 (Caselli et al. 
2002b). The N2D+/N2H+ ratio of 0.05 in L492 is comparable to that in L1498 (0.04), while 
it is lower than in L1544 (0.23) (Crapsi et al. 2005). These results suggest that chemical 
evolutionary stage of L492 would be younger than those of L1498 and L1544, and molecules 
would be less depleted in L492 than in L1498 and L1544. 
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5. Evolutionary Scenario of L492 and Other Cores 

Recently, discussions considering both dynamical and chemical evolution of dense cores 
have been reported (Lee et al. 2003; Aikawa et al. 2005; Crapsi et al. 2005; Shirley et al. 2005). 
For comparison, we summarize the evolutionary stages of L492 and other starless cores on 
the basis of signature of dynamical collapse, degree of depletion, and molecular abundance 
ratios, as shown in Table 6. As already reported, L492 is in the dynamically collapsing phase 
as well as L1498 and L1544 (Ohashi et al. 1999; Lee et al. 2001; Crapsi et al. 2005), while 
L1521B, L1521E, and TMC-1 are not (Lee et al. 1999; Hirota et al. 2002, 2004). Therefore, 
L492 can be regarded as a dynamically evolved core. Crapsi et al. (2005) recently reported 
that a degree of depletion of CO in L492 is comparable to that in L1498 and is lower than 
in L1544. On the other hand, the depletion of CO would not occur significantly in the 
carbon-chain-producing regions such as L1521E (Tafalla & Santiago 2004). Our results on 
the CS and CCS distributions suggest that depletion of these molecules seems to be less 
significant than in L1498 and L1544 (Tafalla et al. 2002, 2004). Therefore, a degree of 
depletion of molecules in L492 is intermediate between chemically young cores (L1521E) and 
evolved cores (L1498 and L1544). The deuterium fractionation ratios of DNC/HN^^C and 
DCO"'"/H^'^CO^ in L492 are lower than those in L1498 and L1544, and are comparable to 
those in carbon-chain-producing regions. This means that L492 is chemically younger than 
L1498 and L1544. Finally, the NH3/CCS ratio in L492 is larger than those in carbon-chain- 
producing regions and is smaller than in L1498 and L1544, suggesting that L492 is more 
evolved than in carbon-chain-producing regions but is younger than L1498 and L1544. 

Considering these results, we conclude that the evolutionary stage of L492 is interme- 
diate between the carbon-chain-producing regions (L1495B, L1521B, L1521E, and TMC-1) 
and evolved prestellar cores (L1498 and L1544). We propose the following scenario for chem- 
ical and dynamical evolution. At the initial state of dense cores before protostellar collapse, 
such as L1495B, L1521B, L1521E, and TMC-1, abundances of carbon-chain molecules are 
extremely high because of high abundance of the neutral carbon atom (Suzuki et al. 1992). 
On the other hand, NH3 and N2H+ are not yet abundant because they are produced by 
slow chemical reactions (Suzuki et al. 1992; Bergin & Langer 1997; Aikawa et al. 2005). In 
addition, the deuterium fractionation in these cores are not enhanced because of the low 
H2D"''/Hj]' ratio at the earliest stage of chemical evolution, and also because of a low degree 
of depletion of CO, CS, CCS, and other molecules (Roberts & Millar 2000; Turner 2001; Hi- 
rota et al. 2001; Saito et al. 2002). As a core evolves, abundances of carbon-chain molecules 
are still high while those of NH3 and N2H"'' increase by the gas-phase reactions (Suzuki et 
al. 1992; Bergin & Langer 1997), and hence, the NH3/CCS ratio increases with time. This 
would be the stage of L492. For L492, depletion of CS and CCS is less significant than 
L1498 and L1544, while the CO molecules start to deplete (Crapsi et al. 2005). Depletion of 
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molecules would be just started in the central part of L492, if the double peak structure of 
the H^^CO"'' ( J=l-0) and CCS (Tg-Gs) maps reflect real abundance variation. At the earliest 
phase of molecular depletion such as in L492, the deuterium fractionation is not significantly 
affected by the depletion. In this phase, protostellar collapse has already started (Lee et al. 
2001). As a core further evolves, the abundances of carbon-chain molecules decrease due to 
gas-phase reactions and freeze-out of molecules onto grains while abundances of NH3 and 
N2H+ increase significantly (Suzuki et al. 1992; Bergin & Langer 1997). This results in high 
NH3/CCS ratio. The L1498 core would be just around this stage. The increase in the deu- 
terium fractionation seems to be different from molecule to molecule. The N2D"'"/N2H+ ratio 
increases more slowly than the DCO"'"/HCO"'" and DNC/HN^^C ratios. At the most evolved 
phase like L1544, depletion of molecule is significant, and the deuterium fractionation ratio 
increases drastically. 

It is also likely that timescales of gas-phase chemical evolution, freeze-out of molecules 
onto grains, and dynamical collapse are different from region to region (e.g. Taurus and 
Aquila Rift). For example, Lee et al. (2003) and Shirley et al. (2005) classified three types 
of prestellar cores; chemically evolved and dynamically young cores (L1512, L1498), chemi- 
cally evolved and dynamically evolved cores (L1544), and chemically young and dynamically 
evolved cores (L1689B). Lee et al. (2003) suggest that L1689B, which is in the Ophiuchus 
region, has the different environment from those of L1512 and L1544 in the Taurus region. 
It should be noted that L492 is also a chemically young and dynamically evolved core out- 
side the Taurus region. In order to understand this regional difference, further systematic 
observational studies and theoretical modeling would be needed. 
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vations. TH thanks to the Inoue Foundation for Science (Research Aid of Inoue Foundation 
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Fig. 1. — Integrated intensity maps of molecular lines in L492. The velocity range of integra- 
tion (in km s~^), the interval of the contours (2cr, in K km s~^) and the lowest contour {3a, 
in K km s'^) are (7.0-9.0, 0.24, 0.36) for ^3CO(J=1-0), (7.0-8.5, 0.070, 0.105) for C^^S{J=1- 
0), (7.2-8.3, 0.100, 0.150) for CCS(43-32), (7.0-8.5, 0.200, 0.300) for HC3N(J=5-4), (7.2-8.6, 
0.100, 0.150) for NH3(1,1), (7.2-8.8, 0.150, 0.225) for C^^O{J=1-0), (7.3-8.2, 0.050, 0.075) for 
C=^^S(J=2-1), (7.2-8.3, 0.060, 0.090) for CCS(76-65), (7.2-8.3, 0.060, 0.090) for HC3N(J=9- 
8), (7.3-8.3, 0.100, 0.150) for HC5N(J=17-16), (7.2-8.3, 0.110, 0.165) for DCO+(J=1-0), 
(7.2-8.5, 0.080, 0.120) for H"CO+(J=1-0), (7.0-8.5, 0.070, 0.105) for DNC(J=l-0), (7.0-8.5, 
0.080, 0.120) for HN^3q(^j^;l-0), and (-0.6-15.2, 0.200, 0.300) for N2H+(J=l-0, total). For 
the NH3(1,1) line, the integrated intensity is the sum of all 5 components (see text and Ta- 
ble 1). For the N2H"^(J=l-0) line, the velocity range of integration is for the Fi,F=2,3-l,2 
component. The beam size (HPBW) is shown in the lower left of each map. The peak 
position of the dust continuum emission (Crapsi et al. 2005) is indicated by a filled square, 
and observed positions are indicated by small crosses. Gray scale maps are superposed to 
the ^^CO(J=1-0) and C^^O(J=1-0) maps to clarify their complex structures. 
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Fig. 2. — Examples of the observed spectral lines. Left: Spectra of sulpher-bearing molecules. 
They are observed toward the (40", 0") position of L492. Middle: Spectra of cyanopolyynes 
and NH3. They are observed toward the (40", 0") position of L492, except for NH3(2,2) and 
HC7N(J=21-20), which are observed at the (0", 0") position. Right: Spectra of deuterated 
molecules and CO isotopes. They are observed toward the (20", 20") position of L492. 
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Fig. 3.— A spectrum of HC3N(J=5-4) at the (40", 0") position of L492. A dashed hne 
represents the best fit model. The expected hyperfine pattern of the HC3N hne is indicated 
in the bottom of the figure. 
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Fig. 4.— A spectrum of N2H+(J=l-0) at the (20", 20") position of L492. A dashed hne 
represents the best fit model. The expected hyperfine pattern of the N2H"'" line is indicated 
in the bottom of the figure. 
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Fig. 5.— A spectrum of NH3(J, K=l,l) at the (0", 0") position of L492. A dashed hne 
represents the best fit model. The expected hyperfine pattern of the NH3 hne are indicated 
in the bottom of the figure. Groups of blended hyperfine components are labeled as adopted 
by Ungerechts et al. (1980). 
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Fig. 6. — Integrated intensity ratio of deuterated molecules as a function of the distance from 
the dust continuum peak. Top: DNC/HN^^C ratio. Bottom: DCO^/H^'^CO^ ratio. Errors 
are estimated from the rms noise of each spectrum {2a). Because these lines are optically 
thin, the integrated intensity ratios approximately represent the deuterium fractionation 
ratios as a function of the distance from the dust continuum peak (see text). 
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Fig. 7. — Channel maps of the CCS(76 — 65) hne in L492. The velocity range is labeled at the 
top right corner in each panel. The interval of the contours and the lowest contour are 0.030 
K km s~^ and 0.045 K km s~^, respectively. The peak LSR velocity of L492 corresponds to 
7.7-7.9 km s~^. The position-velocity diagrams shown in Figure 8 are made along the bold 
solid lines in this panel. 
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Fig. 8. — Position-velocity (p-v) diagrams of the CCS( JAr=76-65) line. The top panel shows 
the p-v diagram along the north-south direction at the offset in right ascension of 20", and 
the bottom panel shows that along the east-west direction at the offset in declination of 0", 
as shown in Figure 7. The interval of the contours and the lowest contour are 0.10 K for 
both panels. 
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Table 1. Gauss fit parameters for the observed lines 





Frequency 




t; 


Vlsr 


Av 


J-rms 


Transition 


(MHz) 


Position 


(K) 


(km s-l) 


(km s-l) 


(K) 


13CO(J=1-0)'' 


110201.353 


(20", 20") 


2.5(12) 


7.92(15) 


1.6(10) 


0.3 


Ci8O(J=l-0)=' 


109782.173 


(20",20") 


1.20(15) 


7.69(5) 


0.62(13) 


0.15 


DCO+(J=1-0)=- 


72039.331 


(20",20") 


0.83(17) 


7.70(5) 


0.49(13) 


0.13 


DNC(J=l-0)'' 


76305.717 


(20",20") 


0.67(17) 


7.67(9) 


0.73(22) 


0.09 


H13CO+(J=1-0)=' 


86754.330 


(20", 20") 


0.83(10) 


7.83(4) 


0.61(10) 


0.08 


HNi3c(J=l-0)'' 


87090.859 


(20", 20") 


0.47(13) 


7.82(10) 


0.72(28) 


0.07 


C34s(J=l-0)=- 


48206.946 


(40",0") 


0.71(17) 


7.82(6) 


0.48(14) 


0.08 


C^^S(J=2-1)^ 


96412.961 


(40",0") 


0.47(9) 


7.74(4) 


0.37(8) 


0.06 


CCS(Jiv=2i-lo) 


22344.030 


(40",0") 


0.69(6) 


7.78(3) 


0.80(8) 


0.03 


CCS(Jiv=43-32)'' 


45379.033 


(40",0") 


1.51(10) 


7.73(2) 


0.48(4) 


0.07 


CCS(Jiv=76-65)'' 


81505.208 


(40",0") 


0.87(20) 


7.71(5) 


0.42(12) 


0.09 


C3S(J=4-3) 


23122.983 


(40",0") 


0.41(6) 


7.77(4) 


0.59(11) 


0.04 


HC3N(J=5-4,F=5-5) 


45488.834 


(40",0") 


0.35(13) 


7.76(3) 


0.81(7) 


0.09 


HC3N(J=5-4,F=6-5/5-4/4-3)'' 


*> 45490.306 


(40",0") 


2.35(16) 


7.76(3) 


0.81(7) 


0.09 


HC3N(J=5-4,F=4-4) 


45492.106 


(40",0") 


0.21(13) 


7.76(3) 


0.81(7) 


0.09 


HC3N(J=9-8)'' 


81881.462 


(40",0") 


1.70(9) 


7.73(1) 


0.46(3) 


0.09 


HC5N(J=9-8) 


23963.901 


(40",0") 


1.00(8) 


7.81(3) 


0.75(7) 


0.03 


HC5N(J=17-16)'' 


45264.720 


(40",0") 


1.23(12) 


7.76(2) 


0.49(6) 


0.06 


HC7N(J=21-20) 


23687.898 


(0",0") 


0.18(4) 


7.65(6) 


0.56(15) 


0.02 


N2H+ ( J, Fi,F=l, 1,0-0,1,1)'^ 


93171.621 


(20", 20") 


0.28(10) 


7.75(2) 


0.44(4) 


0.08 


N2H+(J,Fi,F=l,l,2-0,l,2)=' 


93171.917 


(20", 20") 


0.52(11) 


7.75(2) 


0.44(4) 


0.08 


N2H+ ( J, Fi,F=l, 1,1-0,1,0)=' 


93172.053 


(20",20") 


0.48(11) 


7.75(2) 


0.44(4) 


0.08 


N2H+ ( J, Fi,F=l, 2,2-0,1, 1)=^ 


93173.480 


(20", 20") 


0.59(10) 


7.75(2) 


0.44(4) 


0.08 


N2H+ ( J, Fi,F=l, 2,3-0,1, 2)=^ 


93173.777 


(20", 20") 


0.92(11) 


7.75(2) 


0.44(4) 


0.08 


N2H+(J,Fi,F=l,2,l-0,l,l)=' 


93173.967 


(20", 20") 


0.40(11) 


7.75(2) 


0.44(4) 


0.08 


N2H+(J,Fi,F=l,0,l-0,l,2)=' 


93176.265 


(20", 20") 


0.67(10) 


7.75(2) 


0.44(4) 


0.08 


NH3(J,ii'=l, 1; VH)='.<= 


23692.955 


(0",0") 


0.19(7) 


7.87(4) 


0.86(10) 


0.04 


NH3{J,K=1, 1; H)''''= 


23693.895 


(0",0") 


0.31(7) 


7.87(4) 


0.86(10) 


0.04 


NH3{J,K=1, 1; M)='''= 


23694.496 


(0",0") 


0.74(8) 


7.87(4) 


0.86(10) 


0.04 


NH3(J,ii'=l, 1; L)''''^ 


23695.095 


(0",0") 


0.33(7) 


7.87(4) 


0.86(10) 


0.04 


NH3(J,i^=l, 1; VL)''''= 


23696.037 


(0",0") 


0.24(7) 


7.87(4) 


0.86(10) 


0.04 


NH3(2,2) 


23722.634 


(0",0") 


0.06(2) 


7.71(11) 


0.63(27) 


0.02 



Note. — The numbers in parenthesis represent three times the standard deviation in the Gaussian fit in 
unit of the last significant digits. 

''Maps of the fines is shown in Figure 1. 

''Three hyperfinc components are blended. 

'^Hyperfinc components are blended. Each component is labeled as adopted by Ungerechts et al. (1980) 
and Figure 5. 
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Table 2. Results of LVG calculations 









/.b 






Te, 


n(H2) 


N 




Molecule 


Transition 


S^,^' 


(Debyc) 


Position 


TO 


(K) 


(xlO^cm-^) 


(xlOi-'^cm-^) 


Reference 


DCO+ 


J=l-0 


1.00 


4.07 


(20", 20") 


0.47 


7.4 


O.OO'^ 


0.109 


1,2 


DNC 


J=l-0 


1.00 


3.05 


(20", 20") 


1.9 


4.4 


O.OO'^ 


0.56 


3,4 


Hi3cO+ 


J=l-0 


1.00 


4.07 


(20", 20") 


0.82 


5.9 


0.90^= 


0.137 


1,2 


HN13C 


J=l-0 


1.00 


3.05 


(20", 20") 


2.2 


3.9 


0.90== 


0.44 


3,4 


C^^s 


J=l-0 


1.00 


1.96 


(40",0") 


0.20 


8.6 


0.90 


0.38 


5,6 




J=2-l 


2.00 


1.96 


(40",0") 


0.79 


4.7 








CCS 


JjV =43-32 


3.97 


2.81 


(40",0") 


0.43 


9.0 


0.68 


1.61 


7,8,9 




>/jV =76-65 


6.97 


2.81 


(40",0") 


0.54 


7.1 









''Intrinsic line strength. 

''Dipolc moment. 

'^The H2 density derived from the C'^'^S data is assumed. 



References. — 1: Haese & Woods (1979); 2: Montciro (1985); 3: Blackman ct al. (1976); 4: Green & Thaddeus 
(1974); 5: Winnewisser & Cook (1968); 6: Green & Chapman (1978); 7: Murakami (1990); 8: Yamamoto et al. (1990); 9: 
Wolkovitch et al. (1997) 
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Table 3. Results of LTE calculations 









/.b 






Te, 


N 








Molecule 


Transition 


S^i'^ 


(Debye) 


Position 


TO 


(K) 


(xl0i=^cm-2) 


NOTE 




Reference 


13CO 


J=l-0 


1.00 


0.11 


{20",20") 


4.0 


8.3 


1960 


i3CO/Ci80= 


=7.3 


1 


CISQ 


J=l-0 


1.00 


0.11 


(20", 20") 


0.64 


8.3 


310 


i«co/ci»o= 


=7.3 


1 


CCS 


"fjV =43-32 


3.97 


2.81 


{40",0") 


3.4 


5.0 


5.3 


Tex fixed 




2,3 


CsS 


J=4-3 


4.00 


3.64 


{40",0") 


0.21 


5.5 


0.55 


Tex fixed 




2,3 


HC3N 


J=5-4 


5.00 


3.72 


{40",0") 


8.8 


6.4 


17.2 


hf fit<= 




4 


HC5N 


J=17-16 


17.00 


4.33 


{40",0") 


0.65 


6.5 


4.1 


Tex fixed 




5 


HC7N 


J=21-20 


21.00 


4.82 


{0",0") 


0.03 


10.0 


0.32 


Tex fixed 




6 


N2H+ 


J=l-0 


1.00 


3.40 


(20", 20") 


6.1 


4.9 


0.66 


hf fit<= 




7 


NH3 


(1,1) 


1.50 


1.46 


(0",0") 


2.3 


4.1 


34 


hf fit<: , Trot- 


=10 K 


8 



^Intrinsic line strength. 

''Dipolc moment. 

'^Total optical depth and excitation temperatures are derived from a simultaneous fit to the hyperfine (hf) components. 



References. — 1: Goorvitch (1994); 2: Murakami (1990); 3: Yamamoto et al. (1990); 4: Lafferty & Lovas (1978); 5: Alexander 
et al. (1976); 6: Botschwina et al. (1997); 7: Havenith et al. (1990); 8: Cohen & Poynter (1974) 
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Table 4. Column densities of the selected molecules in unit of 10^^ cm ^ 



Molecule 


L492 


L1495B 


L1521B 


L1521E 


TMC-1 


Taurus 


Others 


C34s 


0.38 




0.56=^ 


1.25^ 


0.73^ 


0.17-2.18=' 


0.08-0. 68=" 


CCS 


5.3'' 


1.44^= 


3.6^ 


2.8<= 


6.6^ 


<0.09-6.6d 


<1. 1-1.51'' 


CgS 


0.55 


0.68<='^ 


1.6<= 


1.46 


1.3'* 


0.42-1. 3^^ 




HC3N 


17.2 


2.1<= 


4.1^ 


2.3s 


17.1<^ 


<0.07-17.1<' 


<0.07-1.35<' 


HC5N 


4.1 


0.52° 


1.2<= 


0.468 


5.6d 


<0.14-5.6d 


<0.17-1.07<' 


HC7N 


0.32 




0.241^ 




1.4^ 


0.08-1.4'' 




NH3 


34 


5.5<= 


12.6<= 


7.3s 


19<' 


<2-107'' 


<3-112<' 


N2H+ 


0.66 


<0.17'= 


<0.19'= 


<0.14g 


0.74' 


0.3-1. 3J 


0.14-2. 7J 


H13CO+ 


0.137 


0.046^= 


0.063^= 


0.106<= 


0.14"^ 


0.038-0.164' 


<0.0110-0.126' 


CiSQ 


310 


190™ 


300™ 


170™ 


330"^ 


60-520' 


30-700' 



Note. — Positions for the cores are as follows: L492: a2000 = 18''15'"46M, <52ooo = -03°46'13" 
and Tables 2 and 3; L1495B: 02000 = 04'»15'"36^5, <52ooo = +28°47'06"; L1521B: 02000 = 
04''24™12».7, 52000 = +26°36'53"; L1521E: 02000 = 04'^29'"16^5, <52000 = +26°13'50"; TMC- 
1; 02000 = 04''41™42».5, ^2000 = +25°41'27". 

^'Hirota et al. (1998). 

''The LTE value is employed here for consistency with other results. 

=Hirota et al. (2001). 

^Suzuki et al. (1992). 

"Hirota et al. (2004). 

Position is (80", 40") offset from the reference position. 

Sffirota et al. (2002). 

''Cernicharo et al. (1986). 

'Hirahara et al. (1995). 

JCaselli et al. (2002a). 

••Pratap et al. (1997). 

'Butner et al. (1995). 

■"Myers et al. (1983). 
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Table 5. Molecular abundance ratios as an indicator of chemical evolution 



Molecule 



L1495B L1521B L1521E TMC-1{CP) L492 



L1498 



L1544 



DNC/HNi^c 


<0.66=' 


0.70'' 


0.66*^ 


1.25= 


1.27 


1.91=" 


a.o'^ 


DCO+/H13CO+ 


im^ 


I.IO'' 


0.63=^ 


0.77^ 


0.80 


2.7° 


3.1-9.2^ 


NH3/CCS 


S.S'^'S 


3.5S''' 


2ga,6 


2.9^ 


6.5 


25'^ 


15h,i 



Note. — Positions for the cores are as follows: L1498: 02000 = 04'"10™51^.5, i52ooo = +25°09'58"; 
L1544: 0:2000 = 05^04"16^6, ^2000 = +25°10'48". 

''Hirota ct al. (2001). 

''Hirota et al. (2002). 

<=Hirota et al. (2003). 

•^Turner (2001). 

''Butncr ct al. (1995). 

fCaselli et al. (1999). 

SHirota ct al. (2004). 

^Suzuki et al. (1992). 

'Position is (-19", 60") offset from the reference position. 



Table 6. Evolutionary scenario of dark cloud cores 



Item 



L1495B 



L1521B L1521E TMC-l(CP) L492 L1498 L1544 



Gas-phase molecular abundance 



DNC/HNi3C=' 
DCO+/Hi3cO+ = 

N2D+/N2H+'' 
NHs/CCS^ 



Y 


Y 


Y 


Y 


Y 


E 


E 


Y 


Y 


Y 


Y 


Y 


E 


E 


? 


7 


Y 


Y 


E 


E 


E 


Y 


Y 


Y 


Y 


I 


E 


E 



Degree of depletion 


Depletion factor of CO'' 

Map structure of CS and CCS= 


? 
Y 


7 

Y 


Y 
Y 


7 
Y 


I 
I? 


I 
E 


E 
E 


Kinematics 


Signature of infall'* 


7 


Y 


Y 


Y 


E 


E 


E 



Note. — Y, I, and E indicates young, intermediate, and evolved, respectively. 

"Table 5 

''Crapsi et al. (2005) 

=Hirota et al. (2002, 2004); Tafalla et al. (2002); Tafalla & Santiago (2004) 

•^Lee et al. (1999, 2001); Crapsi et al. (2005) 



